Anisotropic 2D diffusive expansion of ultra-cold atoms in a disordered potential 
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We study the horizontal expansion of vertically confined ultra-cold atoms in the presence of 
disorder. Vertical confinement allows us to realize a situation with a few coupled harmonic oscillator 
quantum states. The disordered potential is created by an optical speckle at an angle of 30° with 
respect to the horizontal plane, resulting in an effective anisotropy of the correlation lengths of a 
factor of 2 in that plane. We observe diffusion leading to non-Gaussian density profiles. Diffusion 
coefficients, extracted from the experimental results, show anisotropy and strong energy dependence, 
in agreement with numerical calculations. 
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Transport in most materials is determined by the com- 
plex interplay of many ingredients, for instance the struc- 
ture and thermal fluctuations of the substrate Q , the in- 
terparticle interactions, which can induce superconduc- 
tivity Q or metal-insulator transitions [Hj], and disor- 
der [fjj. Disorder is relevant to many condensed-matter 
systems and strongly affects transport via scattering. Its 
primary effect is thus diffusion, an effect underlying the 
Drude theory of conductivity as well as the self- 
consistent theory of Anderson localization Q. Disor- 
der is of special interest in dimension two (2D), which is 
the marginal dimension for return probability in Brow- 
nian motion and for Anderson localization [8j. More- 
over, intriguing effects, which are not fully understood, 
occur in 2D, such as the metal- insulator transitions in 
high-n^ilitv Si MOSFETs d, GaAs^heterostruc- 
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Ultra-cold atomic gases are good candidates to study 
classical or quantum disordered systems (see Refs. [14, 
[l5l | and references therein). They offer unique versatility 
as one can control the amount and type of disorder, the 
interaction strength or the confinement geometry. In ID, 
Anderson localization [l6l 17 1 and interaction- induced 
derealization [lS|j have been observed. In 3D, the com- 
petition between interaction and disorder has been inves- 
tigated in disordered optical lattices 3, 2fj. Diffusion 
was reported for speckle-induced 3D optical molasses in 
the dissipative regime [2l[. So far, less work has been 
devoted to 2D. 

In this letter, we study diffusion of ultra-cold atoms 
in an effectively anisotropic disordered potential without 
dissipation. The geometry is planar as the atoms are con- 
fined vertically to a size of about 1 /xm in a dipole trap 
and horizontally free to move over a millimeter. In the 
presence of disorder, we observe expansion at a reduced 
speed and anisotropic, non-Gaussian atomic density pro- 
files. We show that the dynamics is horizontally diffusive. 
Fitting a diffusive model to the data, we extract the dif- 
fusion coefficients and find that they are anisotropic and 
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FIG. 1: Experimental setup. The ultra-cold atom cloud is 
released in a vertically confining potential created by a blue- 
detuned Hermite-Gauss TEMOl-like beam [H, HU, prepared 
with a 0-7r phase plate (left). The disorder is induced by an 
angled speckle light field created from a diffusive plate (right). 
Finally, the horizontal cloud expansion is detected from the 
top by fluorescence imaging on an Andor Electron Multiplying 
Charge Coupled Device (EM-CCD). 



strongly energy-dependent. Our results are consistent 
with numerical simulations assuming classical dynamics. 

The experimental setup uses a vertically confining po- 
tential and a speckle light field (see Fig.Q]). Both are 
created with 767 nm laser light [22j, blue detuned from 
the resonance at 780 nm for Rb atoms in their ground 
state. They thus induce a repulsive potential. The ver- 
tical confinement is realized between the two lobes of a 
vertically focused Hermite-Gauss TEMOl-like mode, pre- 
pared with a holographic 0-tt phase plate [H, [24| ■ The 
measured vertical trapping frequency is uj/2tt — 680 Hz 
with a 22 fj,m separation between the two intensity max- 
ima, a total power of 150 mW, and an horizontal waist 
radius at 1/e 2 of 1.1mm. 
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The speckle light field is produced by a beam pass- 
ing through a diffusive plate 25, 2(| and focused on the 
atoms. The high numerical aperture (a 75 mm diame- 
ter aperture at a distance of 150 mm from the atoms) 
allows us to achieve an approximately Gaussian corre- 
lation function with transverse correlation length a y — 
0.8 /im (half- width at 1/y/e). The longitudinal correla- 
tion length (in the direction of propagation) is deduced to 
be uiong ~ 9^tm [26|. As the speckle beam is at 30° from 
the horizontal expansion plane, the correlation length a x 
along x is 2a y = 1.6 /im. With a power of 66 mW and 
a Gaussian waist radius of 1.1mm, the average disor- 
dered repulsive potential at the center of the beam is 
V ~ k B x 53(8) nK« (2tt7j) x 1.1(2) kHz (where k B is the 
Boltzmann constant and 2ttH the Planck constant). 

The experiment proceeds as follows. An ultra-cold 
atom sample is produced by an all-optical runaway evap- 
oration in a crossed dipole trap at 1565 nm, as described 
in Ref. 271. The atom cloud is first transferred in 5 ms in 
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a trap combining simultaneously the initial crossed dipole 
trap and the vertically confining beam. The crossed trap 
is then further ramped down in 200 ms in order to re- 
duce the confinement and thus also the temperature to 
k B T « fc B x 200(20) nKw (2nh) x 4.2(4) kHz, slightly 
above the condensation threshold. Finally, the speckle 
field is ramped up in 4 ms, and 1 ms later a thermal cloud 
of N = 1.5 x 10 5 atoms is released in the horizontal plane 
by suddenly turning off the crossed dipole trap. After a 
chosen 2D expansion time, the vertical confinement and 
the speckle potential are switched off. After 0.1ms, the 
atomic column density is measured from the top through 
fluorescence imaging. 

A typical image for an expansion time of 50 ms in the 
disordered potential is presented in Fig. [5^. The corre- 
sponding integrated density along y (resp. x) is plotted in 
Fig.[2p (resp.Hfc). We observe a sharp anisotropic struc- 
ture elongated along x around the initial position, sur- 
rounded by a broader isotropic cloud similar to what 
is observed in the absence of disorder. The sharp 
anisotropic structure corresponds to low energy atoms, 
whose expansion has been slowed down by the disorder, 
whereas the broad cloud corresponds to atoms which ex- 
pand almost ballistically at this time scale. For an ex- 
pansion time of 200 ms (see Fig. 2b and 2c), the contri- 
bution of the ballistic atoms is negligible with respect 
to the lowest energy atoms. The cloud profiles are then 
found to be non-Gaussian with long tails in both direc- 
tions. Similar profiles have been theoretically predicted 
for energy dependent diffusive behavior in the expansion 
of Bose- Einstein condensates [28| . 

We first study the behavior of the peak column den- 
sity n(0, 0,i) as a function of time. It should scale as 
1/t in a diffusive regime, and as 1/t 2 for a ballistic ex- 
pansion. Figure [3] shows a log-log plot of the measured 
peak density as a function of time, with and without dis- 
order. For times below 15 ms, the cloud is smaller than 
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FIG. 2: Atomic column density after planar expansion of an 
ultra-cold gas in an anisotropic speckle potential, a: Image 
after 50 ms of expansion. b,c: Integrated density along the 
two major axes. The plain dots (open squares) correspond to 
50 ms (200 ms) of expansion. 
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FIG. 3: Evolution of the peak column density n(0, 0, t) as a 
function of expansion time. Triangles: with disorder; circles: 
without disorder. The solid lines correspond to fits with al- 
gebraic time dependence between 15 and 200 ms. The fitted 
slopes of the decay are -0.98 with disorder and -1.97 without 
disorder. 



the pixel size and therefore our measurement does not 
reflect n(0,0, t). Between 15 and 200 ms, we observe a 
linear behavior with a slope — 2.0^Q'g without disorder, 
whereas with disorder, we find a linear behavior with a 



slope —1.0 
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The uncertainties come from the disper- 



sion of the slopes found for different data sets taken in 
similar conditions. This measurement is consistent with 
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FIG. 4: a: Diffusion coefficients in the x (blue circles) and 
y (red diamonds) as a function of energy in log-log scale. 
Points are numerical results, lines are fits to power laws 
Di-(E) = D^E/V)^, £ = x,y. b: Energy distribution for 
the experimental parameters. The shaded regions correspond 
to sub-diffusive regimes (see text). 



diffusive expansion of a significant part of the atoms in 
the horizontal plane. After only 15ms, the contribution 
of the ballistically expanding atoms to the density at the 
origin vanishes. 

In order to understand our experimental findings in 
more details, we have performed numerical simulations. 
In the experiment, k^T ~ Qftw, so that a few vertical har- 
monic oscillator states are populated. The vertical size 
of the atomic cloud, Az ~ yj k-gT /moj s» lum (where 
m is the atom mass), is much smaller than cri ong and 
the speckle potential can be considered invariant along 
its propagation axis. Since it makes an angle 9 = 30° 
with respect to the expansion plane, it couples the verti- 
cal quantum states. To account for these features in the 
numerics, we consider the 3D dynamics of classical parti- 
cles in the external potential V(x, y, z) + muj 2 z 2 /2, with 
V(x, y, z) — Vi so (a; sin 9 — z cos 9, y) where V mo (u,v) is a 
2D isotropic speckle potential with correlation length o~ y . 
A classical particle model is reasonable since k-oT/hw ps 6 
and ko~ y ~ 5 where k = ymfeT/fi. 

A characteristic time scale for the dynamics is the 
Boltzmann time tb, i.e. the time after which the mem- 
ory of the direction is lost, which depends on the par- 
ticle energy. For long times (t ^> tb), scattering from 
the angled speckle potential redistributes the kinetic and 
potential energies in 3D, so that the dynamics in the 
horizontal plane is expected to depend on the 3D parti- 
cle energy E. We hence calculate the spatial variances 
(£, 2 {E, t)) as a function of time, where £ = x, y and brack- 
ets indicate averaging over disorder and over initial con- 
ditions corresponding to the energy E. The evolution 
is described by (£ 2 (E,t)) ~ 2D(:{E)t^ E \ We iden- 
tify three regimes characterized by the value of J^(E). 
For E/V < 2, we find a subdiffusive dynamics, i.e. 
75 (E) < 1, for experimentally-relevant time scales. In 
particular, for E/V < 0.52, we find strictly bounded 
trajectories, j^(E) = 0. This is consistent with the 
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FIG. 5: Two-dimensional density distribution after 200 ms of 
expansion. Density profiles through cuts along x (resp. y) at 
different positions along y (resp. x). The upper (red) points 
are cuts through the central peak, whereas the other curves 
downwards corresponds to positions separated by ^36 /im. 
They are artificially offset for clarity. The lines are the re- 
sult from the 2D fit with Eq. (1) convolved by the imaging 
resolution (~ 15 (im) . 



percolation threshold expected for 2D speckle potentials 
|29l [30I ] . For E/V > 2, numerical simulations yield 
a diffusive dynamics, i.e. J^(E) ~ 1. In this regime, 
the diffusion coefficients D^(E) are strongly anisotropic 
and grow algebraically with the particle energy (Fig.0^). 
From a fit to the numerical calculations for our pa- 
rameters, we find D X (E) = 2.\Jv o~l/m(E /V) 2 * and 



D y (E) = 0.65JVa%/m(E/V) 2 - s [3J,|32j. We have also 
done simulations of a classical 2D diffusion in the same 
anisotropic disorder. The various regimes found in the 
3D simulations with vertical confinement are also found 
in 2D simulations at the same values of E/V . In the dif- 
fusive regime, the energy dependence of the diffusion co- 
efficients remains algebraic but with modified constants. 

In the experiment, the observed expansion results from 
the diffusion of atoms with a broad energy distribu- 
tion, N(E). It can be calculated assuming that, be- 
fore abrupt release in the horizontal plane, the gas is 
at thermal equilibrium in the trap plus speckle poten- 
tial. The corresponding energy distribution N(E) cx 

e -B/fe B T (l - e (nfto-E)/V-) ig plotted in pig.gj, 

[33j |. It is fully determined from the experimental pa- 
rameters u), N, V and T. Here, only 6% of the atoms 
are sub-diffusive (E < 2V, shaded regions in Fig. 4b). 
Incorporating their contribution to the diffusive regime 
is thus a small error, and for long expansion time, the 
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column density can be approximated by 

( 



n{x,y,t) 



1 exp 
dE N(E) - 



AD x (E)t 4D y (E)t 
4TT^/D x (E)D y (E) 
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(1) 

We fit Eq. (1) (convolved with our imaging resolution 
~ 15 /im) to the experimental 2D density distribution 
with D X (E) = D x (E/E R ) a and D y (E) = D° y (E/E R ) a , 
D x , Dy, and a as fitting parameters, and the recoil en- 
ergy En — fee x 180 nK as energy scale. As can be seen 
on Fig. El the 2D fit function reproduces the data both 
close to the central peak and in the wings. We find D° = 
3.0(1.5) x 10- 7 m 2 .s _1 , D° = 8.7(4.3) x lO^m^nis 1 , 
and a = 3.3(3). The uncertainties come from the un- 
certainties on the measurements of N, V and T used in 
N{E) and from an observed systematic drift of the results 
as a function of the expansion time 3J|. Experimen- 
tally, the power-law exponent is found to be a — 3.3(3), 
to be compared with 2.8 in the simulation. The ob- 
served ratio of the two diffusion coefficients is 3.45(15) 
when it is 3.7 in the simulation. These slight discrep- 
ancies can be due to the approximations made in or- 
der to derive Eq. (1). At E — E R w fc B T, numer- 
ically, we find D X (E R ) = 1.3(0.6) x Kr 7 m 2 .s -1 and 
Dy(En) = 3.5(1.7) xl0 -8 m 2 .s -1 , where the uncertainties 
come from the uncertainties on the experimental param- 
eters V and a y . The experimental values of the diffusion 
coefficients are thus in quantitative agreement with the 
3D classical simulation considering our relatively large 
uncertainties. 

In conclusion, we have observed and studied 2D diffu- 
sive expansion of ultra-cold atoms in a disordered poten- 
tial. As a result of the effective anisotropy of the speckle 
potential, the diffusion is anisotropic. Fitting a diffusive 
model to our density profiles, we are able to extract the 
diffusion coefficients and find a strong dependence on the 
atom energy, in quantitative agreement with a classical 
simulation for our parameters. 

Understanding the diffusion properties as a function of 
energy (in particular out of the weak scattering regime) 
is a necessary step towards the study of other disorder- 
induced effects in two dimensions, starting with anoma- 
lous subdiffusion [35| and classical trapping under the 
percolation threshold 29|, 3(| ■ By cooling the gas further 
or by reducing the correlation length of the disorder, we 
expect quantum corrections to the diffusion and Ander- 
son localization to show up at the sub- mm length scale of 
the experiment 36]. Moreover, in a 2D degenerate gas, 
the influence of disorder on the Berezinskii-Kosterlitz- 



Thouless transition [37| is especially intriguing. Will the 
vortices be pinned by disorder [H}? 
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